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Abstract. We study a particular broken ray transform on the Euclidean unit 
square and establish injectivity and stability for Cq perturbations of the con- 
stant unit weight. Given an open subset E of the boundary, we measure the 
attenuation of all broken rays starting and ending at E with the standard 
optical reflection rule. Using the analytic microlocal approach of Frigyik, Ste- 
fanov, and Uhlmann for the X-ray transform on generic families of curves, we 
show injectivity via a path unfolding argument under suitable conditions on 
the available broken rays. Then we show that with a suitable decomposition of 
the measurement operator via smooth cutoff functions, the associated normal 
operator is a classical pseudo differential operator of order —1 plus a smoothing 
term with Cg° Schwartz kernel, which loads to the desired result. 

1. Introduction 

In this article, we are interested in the questions of injectivity and stabihty for the 
Euchdean broken ray transform on the unit square in 2-dimensions. Much work has 
aheady been done with regard to the attenuated ray transform ( |Bal04| INov02al 
INov02bl IBS041 IFinSBj IFin041 | Nat01j), and fo r the more general geodesic X-ray 
transform on smooth manifolds f |FSU08[[SU04] ). even in the case of generic families 
of curves. In particular, the main condition for injectivity is, roughly speaking, that 
the union of the conormal bundles over all available curves must span the cotangent 
space of the domain (see |FSU08| ). Such aforementioned results typically make the 
assumption that all rays have no conjugate points. In the Riemannian case, this 
is true if one assumes, for example, that the underlying manifold is simple. But 
of course such a condition trivially holds in the Euclidean case. However, more 
recently Stefanov and Uhlmann considered the case of microlocal invertibility for 
the geodesic ray transform when there are possible isolated conjugate points along 
certain geodesies, which is the setting for fold caustics ( jSU12j ). 

The main tool used to prove injectivity in |FSU08j for the generalized X-ray 
transform, and more generally to recover the singularities of an unknown interior 
distribution /, is analytic microlocal analysis. We employ such ideas heavily in this 
work. However, it should be stated that other results pertaining to the attenuated 
ray transform or geodesic ray transform have utilized different methods. 

The framework of this paper is of a simpler geometric nature, since few results 
have been attained for the problem considered so far. Roughly speaking, if (17, g) is 
a compact Riemannian manifold with appropriate conditions, and E C Oil is open, 
the question is whether the integrals of an unknown function / along all broken 
rays starting and ending on E uniquely determine /. Unfortunately, reflection 
tomography of this nature for non-Euclidean metrics has not yet been treated by 
any previous work. Indeed, the technique of |FSU08| relies on constructing local 
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coordinates about some fixed geodesic 70 in order to show the normal operator has 
nice microlocal properties. But, if one were to reflect from the boundary the family 
of geodesies associated with such coordinates, one no longer has local coordinates 
around the reflected segment of 79. As such, we will focus just on the Euclidean 
case, and in particular we will consider the unit square in M'^. It is the author's 
belief, however, that the results herein are generalizable to at least polyhedral 
domains which tile M", though the computational details may be messy. 

There is, however, some work that is directly related to our problem of interest. 
Eskin proved a uniqueness result for unknown functions in 2-dimensions given 
knowledge of broken rays that reflect off of known convex obstacles in the domain's 
interior ( [Esk04| ). There the weight function is taken to be a smooth electric poten- 
tial V{x). One of the main tools used in the proof is an identity of Mukhometov. 
Furthermore, the recent work of Ilmavirta in |Ilml2) proves an injectivity result 
on the disk for open subsets E of the boundary when the unknown function / is 
uniformly quasianalytic in the angular variable (when written in polar coordinates). 

Finally, we mention some recent work of Salo and Kenig, which provides explicit 
motivation for the consideration of such a broken ray transform in the first place 
f |KS12| ). There the authors consider the recovery of a Schrodinger potential on a 
particular class of non-compact manifolds from partial Cauchy data. This of course 
directly relates to the Calderon problem with partial data. It turns out that the 
question of unique determination of the potential from such partial data would be 
resolved if one can show injectivity of the broken ray transform considered here, 
albeit in a more general context. 

The structure of this paper goes as follows. In S|2] we describe the problem 
and notations and then state the main results. |3] applies the microlocal ideas of 
[FSUOSj together with the path unfolding technique to show injectivity in the case 
of attenuation. In ^J4] we pursue stability of the inverse problem by analyzing 
the normal operator. To this end, §4.11 14.2| and |4.3| all deal with computation of 
the localized normal operator for different choices of neighborhoods of broken rays. 
Finally, we extend the stability estimate and injectivity to perturbations of the 
attenuation cr in Sjsj 
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Academy of Finland while at the University of Jyvaskyla and from the Royal 
Swedish Academy of Science while at the Mittag-Leffler Institute. The author 
would like to thank Mikko Salo for introducing him to this problem and for pro- 
viding helpful comments and feedback. 

2. Statement of Main Results 

We will state the problem setup and definition in the case of a general domain 
n with smooth boundary. However, the main theorems themselves assume further 
that rt is the square. Let be a smooth bounded domain in M" and define r± = 
{{x, 9) S dVl X §"^^ I ± v{x) • > 0} as the set of outgoing and ingoing unit vectors 
on 9ri, respectively. Here v{x) is the outward unit normal vector to d^. We also 
define, for a general open subset E C the set 

(1) ^±{E) := {{x,0) eSn\xe E,±e-iy{x) > 0}. 

Assume that a G C°°{fl x S"^^) and vanishes near dfl. 

We will call any unit speed curve jx,e regular broken ray in Q, if 
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(a) its associated tangent vector field {jx,e{t),ix,e{t)) intersects T-{E) 

(b) it consists of finitely many straight line segments 7x,0,i, 7x,0,2, • • • , 7x,6i,iv 
before hitting E again 

(c) it does not intersect dE, and 

(d) it obeys the geometrical optics reflection law 

(2) %,e,j+i{0) = -%,0,j{Lj) + 2 {%,e,j{Lj) - iy{-fx,e,j{Lj)) ■ %,0,j{Lj)) . 

Here Lj denotes the length of 7x,e,j and the subscripts refer to the fact that 
(7x,e(O),7x,e(O)) = (x,0). 

We use the notation T : r_ ^ r_ to denote the billiard map mapping a vector 
{x,6) G r_ to {x',9') E r_, where x' is the intersection point of the geodesic 'y^.e 
with dQ, and 9' is the reflected direction. We use the maps 7ri,7r2 to denote the 
projections r_ dSl and r_ — )• S"~^, respectively. The broken ray transform 
with respect to the accessible set E, denoted by Icr,E '■ L^{fl) — i^(r_, rfS), is then 
defined by 

Ia,Ef{lxfi) = Ia,Ef{x,0) 

■■= / «'^P -II / (^{T^i°T'^{x,e)+Tn2oT"'{x,e),7r2oT^{x,e))dT] 

• exp ( - / CT(7ri o T^{x, 0) + {t - t)it2 o T^{x, 0), 772 o T^{x, 0)) dr] 

(3) 

• /(tTi O {x, 0) + tTT2 o {x, 0)) dt. 

N{x,9) 

= J2 [wjf]{TTioT\x,e) + tTT2oT^{x,e),7r2oT^{x,0))dt. 

N{x,9) 

= Y / [wjf]{zjix,9)+t9,{x,9),9,ix,0))dt. 
for all regular broken rays "fx,e- Here 

(4) {zjix,9),9j{x,0)) = T^{x + T-{x,9)9,9), {x,9) G Sn 

where Sfl is the unit tangent bundle. Also t±{x,9) = min{t > 0\x + 19 € dfl} 
are the positive and negative travel times from a given unit tangent vector {x, 9) e 
SO,, and N{x^ 9) denotes the number of reflections of the given trajectory before 
intersecting E. The weight functions Wj on SVt are given by 

??) = exp I - ^ / a{zm-j + T0m-j,Om-j) rfr cxp I - / a{y - tt], r]) dr J 

= exp -V/ a{z-m + T0-m,O-m)dT] expi - cr(y - r?7, ??) rfr ) 

The measure = • 0\dS{x)d0 is the standard measure on r_. It is im- 

portant to note that because of our assumption that a vanishes near dQ, this 



4 



M. HUBENTHAL 



ensures that the weight functions Wj are smooth. The reason is that we can cir- 
cumscribe fl with the baU B ~ ^((511)1^) ^-^d then replace Zm-jiyTrj) ~ 
TTi o T"'-^^[y + T^{y,ri)j],ri) with the last point of dB that intersects the line 
Zm-j{y,'n) + i(^m-j{y,v) for i < 0. This point is smoothly dependent on (y,??) G 

Remark 1. Although general polyhedral domains lack obvious reflection rules at 
the corners (some authors ignore such trajectories and some define them in terms of 
limits of regular trajectories), the special symmetrical nature of the square does give 
an unambiguous rule for such reflections. But one still runs into the problem that 
the composite functions Zj(x,0),6j{x,9) are not real analytic for {x,0) associated 
with a corner point. One simple, albeit more restrictive route, is to assume that 
the subset E also includes a small neighborhood of each corner. In that case, all 
broken rays with corner points are just ballistic X-rays. More generally, one can cut 
out the rays with corner points from the measurement operator. However, if one 
assumes the attenuation coefficient to vanish in a neighborhood of the boundary, it 
is still possible to treat rays with corners as we will discuss later. 

It is easy to see that N{x, 9) is piecewise constant, which leads to the unfortunate 
fact that Ig- E in general does not preserve the smoothness of a function. Since we 
are considering a convex domain (flat) in this work, the discontinuities of N(x^ 9) 
occur only for {x,9) E r_(i?) such that 'jx,e has a reflection occurring at some 
point in dE. This is the reason they are not classified as regular broken rays. But 
fortunately, this problem can be averted by omitting such rays from the data with 
a smooth cutoff. 

To eliminate the discontinuities introduced by N{x, 9) we will use a collection of 
smooth cutoff functions ak on r_(i?) so as to remove the bad parts of Ia,E- For 
example, for each point (a;o,0o) G r_(i?) where N{x,9) is continuous (i.e. jxo.Oo 
and nearby rays do not touch corner points), we may multiply Ia,E by a smooth 
function a{x,9) = 1 near (xo,9q) and such that N{x,9)\supp{a) is constant. This 
localized operator will have a constant number of reflections and thus behaves very 
much like the usual generalized ray transform, which makes it possible to apply 
the microlocal analytic techniques of _FSU08 . Of course, additional care is needed 
to deal with broken rays that contain corner points, which we discuss later. The 
central idea is that for fixed (x,^) G T*r2 \ 0, any broken ray 7 passing through x 
with direction normal to ^ will allow for the detection of {x,£,) in WF^(/). From 
hereon we use the notation WF^(/) to denote the analytic wavefront set of /. 

The main idea we use in order to apply microlocal analytic techniques to the 
operator is a well-known method of billiards called path unfolding. Essentially, 
we can unfold any broken ray in Q into a straight line on a different domain 
defined by reflecting ft across each of its edges repeatedly to obtain a tiling of M^. 
Similarly, we can extend / to a function / on by reflection across dfl. Finally, if 
K C ft, then we construct K C by concatenating reflected copies of K. We will 
use indices (^1,^2) € to denote how many horizontal and vertical reflections a 
given broken ray undergoes on the unfolded domain before returning to the set 
E. We call this the reflection signature of a path. For example, (—1,2) indicates 
1 reflection to the left and 2 reflections upwards. The following maps and their 
respective inverses will be useful when translating between vectors on Tfi and TO. 
Given a reflection signature (^1,^2); we can compute the reflected version of a point 
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Figure 1. 2-reflected unfolded paths from x to y. 



w S ri by 

(5) 

h + — + (-i)'^w'i,^2 + — + (-l)'^^2 

(6) 

w = R-]^iw) = (^i-iy^+'h + ^^^^ + i-iy^wu i-iy^+'h + ^^^^ + (-i)'^^2 ) . 

Similarly, given an angle 0, Si-^j2{9) is the angle after an (li^h) reflection. 

(7) r,^5z,,,(0) = (-l)'^+'^0+|[l-(-iy^] 

(8) = Si;],irj) = {-ly^ f (-l)S, + 



Such formulas can be verified with straightforward geometry. 

However, applying such unfolding requires significant symmetry conditions on 
the attenuation a. First we define an extension of a on fl x §^ by 

(9) ^ix,0) ^ c7{Ri[]^{x),Sr^\^{9)) iov xe[h,h + l] X [?2,;2 + l], (Zi,?2) eZ2. 

We then require that w{x,9) := eyi-p j^'^^"'^"^ '<7{x — t9 ,9) dr^ is real analytic 

for all {x,9) S x such that the unfolded broken ray through {x,9) is regular. 
We've used the notation L{jx,-e) to mean the length of the half broken ray with 
initial condition (x, ~9) before hitting some copy of E in f2. Obviously, a = 
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works. But such a condition implies that a has symmetries somewhat hke to that 
of a doubly periodic real analytic function, which would necessarily be constant 
(the Weierstrass p- function is doubly periodic but has poles of course). Having 
(7 = near dVL also allows for w to be real analytic since then the singularities of 
t{x^O) near the corners are masked. Thus, for the purposes of this paper, we will 
start with the case that the weight function is identically 1 (equivalently cr = 0), 
and then prove stability for perturbations of a that remain in a neighborhood 
of dVl. 

Note that one can also write Ia,Ef in the form 

(10) i^Ef{x,e)^ w{x + te,e)f{x + te)dt 



Jo 

though we will not use it extensively in the subsequent discussion. 

For a given Nmax G we define the visible set by 
(11) 

Ai :— {x £ £ E>^ , {x, ^) e A^*7 for some regular broken ray 7 with N{j) < N„ 

We can now state the following injectivity result for the 2-dimensional square with 
Euclidean broken rays: 

Theorem 1. Let E C dVl he open and take ct = 0. Let Nmax & let A4 be the 
corresponding visible set, and let K ^ M. Then Iq e is injective on L^{K). 

One can also obtain a stability result analogous to that of jFSUOSj . To do 
so, we must first parametrize a family of regular broken rays, having an upper 
bound on the number of reflections, and whose conormal bundles cover T*n. Fix 
a positive Nmax S N. For each regular broken ray Jxafio with N{xo,Oa) < N^ax, 
let (Zi, ^2) be it's reflection signature on VL. If 7x0,90 1^^^ corner points, then find 
a maximal connected neighborhood U C T-{E) containing [xq^Oq) such that all 
broken rays with initial data in U are regular, have no corner points, and have the 
same reflection signature {li,l2) (See Figure[2|. In particular, U itself contains no 
corner points. We then choose a smooth cutoff function ak on r_ that is equal to 
1 on a compact subset of U and vanishes outside of U. 

If xq is a corner point, then we consider all directed lines C ~ {{x + t9, 9)\t ^M.} 
for (x, 9) in some neighborhood of {x^, 9q) such that N{x, 9) is constant. Then take 
the intersection of C with the ball -B=((^,^),^)of radius ^ whose boundary 

dB circumscribes $7. We then define a smooth function a on dB x §^ supported in 
dB X §^ n £. We may then extend a to x §^ to be constant on all directed lines 
[x + t9, 9). Its restriction to r_(90) is the cutoff we seek. We will treat the case 



that "fxofio contains an additional corner point in ^4.3 



Let a = 'Y^kZi be a sum of smooth cutoff functions on r_(_E) in the sense of 
the previous paragraph, and consider the regularized operator 



E- 



(12) Ia,E,a ^ Ctkla, 

k=l 

Then 

K 

(13) N^^E,a '■= K,E,a^<y.E,a = ^ Ia,E^ki<^k2lcr,I 

ki = l,k2 = l 




Figure 2. A typical beam of unfolded broken rays which follow 
the same reflection sequence before hitting E. The beam termi- 
nates in the part of f2 with reflection signature (/14, V4) = (2, 2). 



The transpose is taken with respect to the standard measure \v{x) ■ 9\dS{x) dO 
on r_(i?). Note that we may as well assume that the ak functions have disjoint 
support, since it won't change the kind of terms present in the normal operator (if 
two such cutoff functions a^^ and have overlapping support, then their product 
is also smooth with the supported broken rays localized just as in the previous 
construction). 

We also define the microlocally visible set by 

(14) M' :== {(z, S)<^T*VL\ (z, £ N*-f^^g for some (x, 6) with a{x, 6) > 0}. 

Remark 2. One reason we need to impose an upper bound on the number of 
reflections is that as N increases, the open level sets of N{X, 0) shrink to zero. 
Similarly, the neighborhoods in r_ on which we have constant reflection signature 
{li,l2) become infinitesimal as + IZ2I — >■ 00. Thus a will not be smooth unless 
we apply some threshold. 

We have the following stability result analogous to that in [FSUOSj . 

Theorem 2. (a) Fix a that is in a neighborhood of dfl, and choose a S 
C°°(r_) as above with reflection threshold N„iax- Fix a set K <s M com- 
pactly contained in the visible set M.. If la, E, a is injective on L^(K), then 

(15) ^II/IIl^(K) < \\Na.E,af\\HHn) < C\\f\\mK)- 

(b) Let a*^ be as above related to some fixed ctq. Assume that Iao.E,a° is in- 



jective on L (K). Then estimate (15) remains true for (it, a) in a small 



neighborhood of (ctq, oP) such that cr = near dft, and with a uniform 
constant C > 0. 
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3. Injectivity of Io,e for the Square 

For the 2-dimensional square, we will be able to establish injectivity of the broken 
ray transform for the unknown function / supported on some subset of the domain 
that depends on E. We will be able to detect covectors from N*jx,e in WF^(/) for 
(most) {x,9) e r_(i?). The idea for the proof of the theorem is heavily borrowed 
from the microanalytic techniques utilized in [FSUOSj . We will need to assume that 
/ is compactly supported within Q, which in particular will allow us to treat the 
reflected ray as a smooth regular curve in the sense of |FSU08) since the integrand 
vanishes near the points of reflection. By using the path unfolding approach, we 
can relatively easily obtain the following microfocal result. 

Proposition 1. Let f € V'{K) for K n. Suppose h,Ef{x,0) = for all {x,0) 
in a small neighborhood of {xo,9o) € r_(_E) where ^xo.9o '■^ ^ regular broken ray. 
T/ien WF^(/)n7V*7,„,eo =0. 

Proof. We need to construct coordinates near the broken ray 79. To do this, we let 
{xq.ii • • ■ , a^cAf} be the ordered reflection points on dfl for 70, where N — N{xq, 0q). 
We then construct a new domain ft consisting of TV reflected copies of fl and glued 
together with the original along the edges that contain each xqj, I < j < N. 
Specifically, let Pm for 1 < m < 4 be the lines in which bound the sides of 
H, in some order. Let {mjJjL^ be the indices corresponding to which edge of dH. 

that xo,j belongs. We iteratively construct Q by starting with xq^i and reflecting 
il across edge P„jj . Then glue the two copies of fi together along . Continuing 
in this way, we obtain a space consisting of + 1 copies of fl, each reflected along 
a specific direction. If 70 contains a corner point, we can add additional refiected 
copies of n as necessary to accommodate nearby broken rays. We also obtain a 
distribution / G 2?' ($7) and an unfolded version of 70, denoted by 70, using the same 
process. Note that by reversing the refiection sequence, any covector (z, ^) G T*ri 
corresponds to a unique covector (z,^) G T*f2. 

Since the edges of are flat, it follows that 70 is a straight line. The openness 
of E and the assumptions on {xo,9o) ensure that for {x,6) in a neighborhood 
U C T^{E) of {xo,9o), the unfolded paths jx,e fie in the same extended space f2 
and avoid dE. Now choose a point po on the line xq + t9Q for t < so that po ^ fl. 
Define x — po + tO. Then are local coordinates near any point of 11 n 72:0.60 • 
We can assume without loss of generality that 6q corresponds to the vector (0, 1). 
Writing x = — {9^,t), we have the coordinates 

U = {x\\x^\ < e, r < t < 1+} C Q. 

Since / is compactly supported inside of fl, we may shrink e as necessary so that 
we can take /~, l~^ to be constant. Now as in the proof of Proposition 1 in |FSU08| . 
we let (zoiCo) G N*jQ. Then consider the corresponding refiected covector via 
path unfolding given by (z*^,^o) G ^*7o- The argument of |FSU08| shows that 
(z*^, ^0) ^ WF^(/). By undoing the reflection process again, we have that (zq, ^0) ^ 
WF^(/). ' □ 



Proof of Theorem^ Suppose la.sf = 0. By Theorem [T] we have that / is analytic 
on K. Since supp(/) C K , we have that / can be extended to an entire function 
and hence must be identically zero. □ 
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Remark 3. In two dimensions, Theorem [T] is not so remarkable since in this case 
every unit vector has a unique unit normal vector up to sign. For example, in the 
case of the square Theorem[T]gives injectivity on the whole domain if E contains two 
adjacent (closed) edges. However, one must be careful since it fails for E consisting 
of two opposite edges. 



4. Stability 

In order to prove the stability result for e with a a perturbation of that 
vanishes near dfl, we will first compute its normal operator. Suppose there is an 
upper limit Nmax on the number of reflections we want to allow in the data, which 
is implicit in the choice of cutoff functions ak{x,9). We construct the modified 
operator 

(16) i„^E,af{x,e):=Y,Mx,o)i^,Ef{x,e), {x,e)er^{E). 

k = l 

We will always suppose that K A4 is fixed, where M C fl is the visible set. We 
will prove the following proposition regarding the structure of Na E,a'- 

Proposition 2. N„^E,a = No- 

,E ,a., ballistic ,E,a,reflect where N^ e, a, ballistic 0, 

classical pseudo differential operator of order —1, elliptic on M, and Na-^E,a,re fleet 
is an operator with Co°(fi x fl) Schwartz kernel. Thus there exists a classical pseudo 
differential operator Q in fl of order 1 such that 

(17) QN^,E,af = f + QN,,E.a,refleetf + Si 

for any f € V(K), and an operator Si with C(^{Q x J7) Schwartz kernel. 

4.1. Computing the Normal Operator . Recall that given j E Z we have 
defined the variables Zj, 9j depending on a; G fi, € S"~^ by 

(18) ■.^T\x + T^{x,0)9, 9) 

Notice that {zo,9o) = {x + t_{x,9)9,9). 

Choose a cutoff function a{x,9) on r_ such that al^^E has a constant number 
of reflections > on its support. Later we will consider broken rays with corner 
points, but then the number of reflections is not constant, so it must be treated 
carefully. The adjoint to aIa,E can be computed as follows 

(19) ^ a{x,9)[UEf]{x,9)g{x,9)dE 

N 

a{x,9)g{x,9)Y^ / Wj{zj + t9j,9j)f{zj + t9j)\v{x) ■ 9\dt dS{x) d9 

For each j we make the change of variables (x,9,t) i~> (y,'?) where {y,r]) := 
{zj{x,9) + t9j{x,9),9j{x,9)). Thus 

{x,9) = {z^j{y,rj),9^j{y,rj)) 

(20) t = -r_(y,77). 

Intuitively, in moving from {y^rj) to {z_j{y,r]),9-j{y,r])), we first project from the 
point y to the boundary along the direction —r] and then trace backwards along 
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the trajectory for the previous j reflections. We compute 

(^Zra{x, 9) + s6l„(x, 6*), 6'„(a;, 0)) 

-3 (y. '7), S-j (y, 7])) + sdrn [z^j (y, rf) , O^j (y, ry)) , 6I„ (z-j (y , r;), 6l_j (y, 77) 



Given a fixed number of reflections j, it is clear that each point {y,rj) G fix S"^-'^ 
corresponds to a unique point {x,9) S r_. Thus, the previously described change 
of variables is a diffeomorphism. The resulting integral is over the domain x S"""'^. 
We have 



N 

E 

3=0 

where 



o^i^-j (y^ v) , S-j (y, v))9{z-j {y, v) , ^-j (y, (y, vi)f{y)Jj {y, v) dy dv, 



Jjiy,v) = 



dxdddt 



dy drj 



Since the billiard map T preserves the symplectic form • 0\dx A d6 on r_, it 

follows that Jj(y, ry) = 1. We stress that this is even true for domains with piecewise 
smooth boundary as is the case for the square, by measure theoretic arguments 
(see |Tab| §1.7). One can verify that the Jacobian factors above cancel by first 
making the change of variables {x,9) T^{x,9) in (19 1, using the invariance of 
the symplectic form \i'{x) ■ 9\dx A d9 under T. Then change to (y, rj) coordinates on 
X S"~^ using the same approach as in (Theorem 1 (b), |SU08) ). Thus we have 



(21) 



N 

(al^.ETgix) = V / [ag]{z^j{x,9),9^,ix,9))wj{x,9)d9 



We are now ready to compute the normal operator /* E^(y,E- By definition 

^a,E,ada,E,af{x) 
N 

= V / \a{z_,,{x,9),9_,,{x,9))\'l„,Ef{z-jAx,9),9_,,{x,9))w,,{x,9)d9 



\a{z_j^ {x, 9),9_j^ {x, 6l))p wj^ (x, 



N N 

EE 

• {zj2-3i {x; 0) + t9j^^j^ (x, 9),9j^^j^ {x, 9)) f {zj^-j, {x, 9) + t9j,_^j^ {x, 9)) dt d9 

N N 



E E 

Jl=0 J2=0j2#il 



(x, 9) Wj^ (zj2-ji [x, 9) + t9j^-j, [x, 9), 9j^..j, [x, 9)) 



■ \a(z_j, {x, 9),9.j, {x, 0))\'' f (z,,_,, (x, 9) + t0,,_,, {x, 9)) dt d9 
+ E / / Wj^{x + t9,9)wj^{x,9)\a{z^j^(x,9),9^j^{x,9))\'^ } {x + T^{x,9)9 + t9) dtd9 
(22) 

= : Nc,^E,a,reflectf{x) + N^^E ,a,bal-listicf {x) ■ 
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It follows that the full normal operator N„^E,cn corresponding to the partition by 
smooth cutoff functions a^, is given by the sum 

No. Nk Nk . „ 

We pick off the ballistic terms where ji = j2 to decompose Na-.s.af as a sum 

.E ,a ^ballistic ,E ,a.,reflect 

with 

,E,a,ballisticf 

= EE/ / \akiz-j,9^,)\^w,{x,e)w,{x + iT.{x,e)+t)9,9) 
■ f{x + {T^{x,9)+t)9) dtd9 



EE/ / \ak{z-j,9.j)\^Wjix,9)wj{x + t9,e)f{x + te)dtd9 
fe=ij=o-^S"-'-^R 

EE/ / [\c^kiz-,ixr)e-,ix,-)mfw,ixr)w,ix + te,-)]^^^J9) 



k = l j=0 

■ fix + t9) dt d9 

and 

Nc Nk r r 

N^,E,a,reflectf{x) =^Y^ ^ / / | Q!fc (^-ji , ) P 

fc=iii=0j2=0j2#ii ' ■^^+ 
• w'ji (x, 9) Wj^ {zj^-j, + t9j^^j^ , Oj^^j^ ) f {zj^-j, + t9j.,-j, ) dt d9. 

The notation [g{x^ ')]even (^) even part of g in the variable 9, defined by 

By Lemma 2 of [FSU08) N^^ E,a,baiiistic is a classical pseudo differential operator 
of order —1 with principal symbol 

K Wfc 

(23) ao(a:,0 = 2^EE / |afe(z„„ 0_,)n«;, (x, d0 

The bound on the number of reflections Nmax ensures that the integral kernel has 
a positive lower bound on the set of all covectors (x, ^) such that it is non vanishing 
at some 9 € S""^ normal to ^. Clearly, Na^E,a,baiUstic is also elliptic on the set M.' . 

4.2. Simplifying Na-^E,a,refiect for the Square . For the reflected part of the 
normal operator, our goal will be to make a change of variables so that /(zj^-ji + 
t^j2-ji) becomes f{y), with the hope that what we obtain is either a weakly singular 
integral operator (i.e. smoothing of order —1) or an operator with Schwartz 
kernel. It turns out that the latter will hold. This technique is very much inspired 
by the kind of substitution used in |SU08[ IHubllj . and also in |FSU08j for the 
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non-Euclidean case. For now, we will still be working with the part of the operator 
that avoids broken rays with corner points. 

Let {x,6) g ri X S^, suppose we have some smooth cutoff a on r_(£'). By 
construction, all broken rays originating from the support of a have the same 
sequence of reflection signatures. Let {{hi, vi), . . . , {h^, vn)} be the ordered set of 



all such pairs (see Figure 2 L Let z 



{x, 9) (i.e. z is the reflection point on 



do, of the broken ray starting at (z^j^ {x, 9),9-j^ (x, 9))). Then (hj^ — hj^ , vj^ — Vj-^) 
is the reflection signature from x to z, and so we consider 

as the associated unfolded coordinates of z. Up to a rotation of fi, there are 
essentially two possibilities depending on whether z lies on a top or right edge of 
a reflected copy of the square. We will use 9 to denote an angle value in radians 
and 9 = (cos 9, sin 9) to denote its associated direction vector. Furthermore, if 
X e E^, we define — {—x'^,x^). For simplicity of notation, we let (lijh) = 

We have one of the following equations for z as a function of the starting point 
X with direction 9: 



(24) 
(25) 



z = (z^^a) = {x^ + {h - x^) cot 9,12) 
z = (/i,i^) = {h,x'^ + {h - a;^)tan6'). 



In the integral expression of the normal operator, we first make the substitution 
V — 'S'ii. (2(^)1 which is the refiected angle after applying the billiard map j2 — 



ji times to {x + T_{x,t 



Note that \li \ + {hi = 32 — ji- Obviously, the 



'32- J 



Ax,9) 



Jacobian determinant of this transformation is 1. We also rewrite 
TTi o T^^~^^{x + T-{x,9)9,9) as the function z{x,Tf). In particular, we may have 
either 
(26) 

z(x, ry) - R-\ [x' + {h - x') cot (^"^ J77)) , ^2) 
(-l)'i+i/i + 



or 
(27) 



zix, = R-]i^ [h, x^ + ih - x') tan [S-],^ (rj)) ) 



1 -(-!)' 



-, (-l)'^+i?2 + 



l-(-l)' 



+ (-1)'^ [x' + ih-x')t^n{SrXiv))] 



Thus we have an integrand involving f(z{x, r])+tfj), where rj — (cos 77, sin 77). The 
next step is to make the substitution y = z{x,ri) + trj (see Figure [s]). We compute 
it's Jacobian determinant to be 

= \{-iy^{l2~'x'^)csc7] + t\ = \{l2-x^)csc9 + t\ = \y-z\ ^ — ^ + 1 



dy 
drjdt 



(28) 



\y-z\ 



— 9 9 

y — X 
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Figure 3. The basic setup for how we define coordinates in the 
integral formula for Na-,E,af- Here j2 = 4 and ji ~ 2. 



if z{x,ni) is given by (26). Similarly, we obtain 
dy 



dr]dt 



y 



when z{x, rj) is given by (27). Here we've used the fact that sec 77 = (—1)'^ seed and 
csc?7 = (—1)'^ csc9. In order to fully carry out the substitution in the integral, we 
need to write z as a function of x and y. We will need the auxiliary equation 

(29) {y-I)-{I~x)^^0, 

which is a consequence of the reflection law. In the case where is undetermined, 
we have z^ = h and then 

(30) t:^ _ & -x^Wi-V^) 



~2 ~2 

z -y 



Similarly, if is undetermined then 



(31) 



{y^-x^m-y") 



y2 



In the given case z^ = h, the Jacobian determinant 



dy 



drjdt 



simplifies to 

1 



y 



X 
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The same holds when — l^. Conveniently, the |y — x\ is strictly positive for 
all y, a; e supp(/) so long as or li is nonzero (i.e. ji — Ji 7^ ). Thus the 
reciprocal factor — resulting from the change of variables is smooth and bounded 

as a function of y and x. Using again the fact that cos 77 = (—1)'^ cos 6', sin 77 = 
(— 1)'^ sin 61 we have that 

^= ((-1)'^ cosf?, (-1)'^ sin0) = ^^((-l)'i(yi - xi), (-l)'ny' ~ 

Altogether, we obtain that N^^E,a.refiectfix) consists of a sum of terms of the form 

(32) / \ak{T-'''{y + T_{y,^)^,^))\^Wj,[x,— : j (y,^) f (y) — r 

V \ \y-x\J \y-x\ 

Evidently, J— is a smooth non- vanishing function of x and y since |y — xj > 
d(supp(/),917) > 0. 

Note that one must be very careful in the above formulation in order to show 



that the integrand in (32 1 is smooth where it is supported. Consider the case if ak 
is non vanishing for a broken ray which starts at a corner point, but still none of the 
broken rays in the support of a contain another corner point after some positive 



number of reflections. Within the factor of (32) given by 

I ttfe (T"-'' (y + T_ (y , 77)77, ) p , 

we see that T~^^ {y+T- {y, r])r], rj) fails to be smooth for (y, 77) that are pulled back to 
this corner point. However, as previously discussed, the construction of ensures 
that it is smooth on x and recall that i3 = i?((i,i),^)is the ball that 
circumscribes fi. By a possible rotation of fi, we may as well assume that hj,Vj 
are all nonnegative. Clearly each sequence {hj}^ {''^j} is also nondecreasing in this 
case. Next define the function qj{y, rj) for I < j < N as the final intersection point 

of the line Rh-^Vj{y) + ^'S';7/t,j '^i*^ for t <Q. We also define qo{y,r]) — y + 
T-.dsiy, 77)^ where T^.dsiy, v) is the negative distance to dB from 77 in the direction 
—77, defined for (77,77) G S x We can then replace ak{T^^^{y + T^{y,r])rj,r])) 
with 



(33) akiqj^{y,v), S^^^^^^^ir])) 



Clearly, (33) depends smoothly on 27,77. 



We observe that the definitions of Wj^ , Wj^ also involve t„ , which is not smooth 
on X However, since a is compactly supported inside of ft on some smooth 
subdomain, we can modify r_ inside of Wj-^ , Wjr^ as necessary to be smooth without 
changing the integrals. In particular, we can use the functions qj to define the 
weight functions instead of the billiard map T and boundary distance function r_ . 



Therefore, (32) is a smoothing integral operator. 



4.3. Incorporating Corners. While the corners of 9r2 have no well-defined nor- 
mal vector, we can make sense of a broken ray that touches a corner point by 
realizing it as a limit of broken rays that avoid the corner. This is described in 
detail in the 2-dimensional case in the work of Eskin in |Esk04j . Essentially, in 2- 
dimensions there are two reflection directions at a corner point that can be obtained 
as a limit of broken rays starting at some fixed point x G and with reflection point 
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approaching the given corner. In higher dimensions, one can obtain a cone of pos- 
sible reflection directions, which is thus a hypersurface. However, for the square 
(or even the n-cube), these limiting rays all conveniently coincide. To see this, let 
X e il and let c be a corner point. Define 9c as the angle between c — x and the 
horizontal axis. For 9 close to 9c, the broken ray from x in the direction 9 will 
intersect twice near c and return at an angle tt + with the positive horizontal 
direction (i.e. it will return in the direction —9). Therefore, the limit of such rays 
shows that "fx,e^ unambiguously reflects back towards x. This is also consistent 
with what one expects from path unfolding. 



Lemma 1. A regular broken ray starting from E contains at most two comer 
points, and if it has two, then the first one is its starting point. 



Proof. Suppose we have a broken ray ^xofio with two distinct corner points. Since 
reflections from corners return in the opposite direction from the direction of inci- 
dence, any broken ray between them will repeat periodically. Thus either one of the 
corner points lies in E or there is a point along the broken ray between them that 
lies in E. But this implies the broken ray would have terminated before hitting 
both corners. □ 



We will use c ^ Xq to refer to a corner in the broken ray if it has one, and we 
will allow for the possibility that xq itself it a corner point. However, as previously 
discussed, xq being a corner point does not present any problem with regard to 
the smoothness of the integrands involved in the expansion of the localized normal 



operator (see the discussion following (32) ). As before, let K CI he the set on 
which we assume / has support. 

So given our broken ray "fxafio > assume it passes through a corner point c ^ xq 
at the m*^ reflection, where 1 < to < — 1, and iV — 1 is the total number of 
reflections of the curve. Choose a neighborhood U consisting of {x, 9) € r_{E) that 
are sufficiently close to {xo,9o) in the following sense. At the m}^ reflection, most 
nearby broken rays to jxo,9o ^il^ undergo a reflection along a very short segment 
near the corner on which / vanishes and cr = 0. In particular, such segments do not 
affect the weights Wj and so we will be able to discard them. Thus, we choose U 
small enough so that at the wth reflection, all nearby broken rays to Jxq.Bq either 
pass through the corner c or hit close enough to the corner so as to reflect into a 
short path that contributes nothing to the integral. Note that for such broken rays 
from U passing very near to the corner c (but not actually touching it), after the 
(to -|- l)*'^ reflection they will be oriented in a direction opposite their respective 
directions during the (to — 1)*'^ segment. We also choose U small enough so as not 
to introduce any new corner points among the included broken rays. 

Given such an open set U, we let lu{x,9) be a function on x §^ which is 
constant on lines and satisfies 
(34) 

1 if {l'x.0(t),'^x,e{t) passes through U and intersects the corner c, 
otherwise 



luix,9) 
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Also let a{x,9) be a smooth cutoff function in x §^ which vanishes outside the 
set of lines associated with U. We write the localized broken ray transform as 



ia,E,af{x,e) = J2[ (1 - iuix,e))a{x,e)[wjf]izj{x,e) + tej{x,e),ej{x,e)) 

N-l . 

(35) +^2 ^u{x, 0)a{x, e)[wjf]{zj{x, 9) + t9j{x, 9), 9j{x, 9)) dt. 



3=0 

We then have 



n,E,M = Yl I [{l-lu)ag]{z.i{x,9),9.^{9))w^{x,9)d9 

JV-l . 

(36) [Iuag]{z_{x,9),9_j{x,9))wj{x,9)d9. 



j=0 

Finally, the normal operator is given by 

N N 



n,E,J.,E,af{x)=J2T. / / [{l-lu)\af]{z-jAx,9),9.jAx,e))wj,{x,9) 
• [Wjj]izj^-n{x,9) + t9j^-j^{x,9),9j^-j^{x,9))dtd9 

N-lN-l „ „ 

+ EE/ / [iu\anz-jAx,0),e-jAx,0))wjAx,0) 



ii=oj2=o 



(37) • [wjj]iz,,_,,{x,9) + t9,,_j,{x,9),9,,_,,{x,9))dtd9. 

Here we've used the fact that (1 — lu)'^ = 1 — Ic/- Then break off the part with 
j2 = ji to obtain 

Na,E,aMiiisticfix) = y2 I I [{l-^u)\aWz-j{x,9),9-j{x,9))wj{x,9) 
~Tq Js^ Jm+ 

• [wj f] ix + T-{x,9)9 + 19, 9) dt d9 

N-l 

+ E/ / [iu\an{z-jix,9),9_^ix,9))wj{x,9) 
•'^^ •'^'^ 

■ [wjf]{x + T-{x,9)9 + t9,9)dtd9.. 

Now we have to be a bit careful in order to combine both summands. In the first 
summand, whenever {z-m{x, 9), 9-m{x, 9)) is in the support of (1 — l[/)a^, we know 
that X + T_(.x, 9)0 + 10 is a line segment very near the corner c where / vanishes. 
Thus we can throw out the term involving j = m. Now for the terms with j < m, 
we may immediately combine the summands into one, thus eliminating the presence 
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of 1(7. That is, we have 

m— 1 



^ E ^a, ballistic 



j=o •'^^ 

• [wjf]{x + T^{x,e)e + te,9)dtd0 

N 

^ / / [{l-lu)\a\']{z-,Ax,0),O^,A^,0))w,,{x,e) 

R+ 



ji=m+l 



Ki /] ix + T^{x,9)e + to, 9) dt d9 

[lu\a\^] {z-j2 {x, 0), 0-32 (2;, &))wj2 {x, 9) 
(38) • [wjJ]{x + T^{x,9)9 + t9,9)dtd9. 



N-l 

E 

j2=m 



For the second and third terms of (38 1, we use the fact that for each j2 with 
m < j2 < N—1, x+(r_ {x, 9)+t)9 in the third term has the same reflection signature 
(li, I2) as in the second term, but with ji = j2 + 1. This has to do with the supports 
of the corresponding cutoff functions. Fix ji,j2 such that m < j2 = ji — 1, and 
let {lijh) be the reflection signature of {x,0) relative to {z_j^{x,9),9_j^{x,9)) in 
the second term or {z-j^{x,9),9-j^{x,9)) in the third. The key idea is that both 

reflection signatures will be the same. Recall the circumscribed disk 5), 
Similar to the case with no corner point, we consider the sequence of reflection 
signatures {{hi,vi), {h2,V2), ■ ■ ■ , {hN-i,VN-i)} in the unfolded domain where all 
of the unfolded broken rays from U touch and over which / has nonzero integral. 
Thus we've cut out the tiny "bow-tie" region near the corner c, which is why there 
are only iV — 1 such pairs. Also recall for < m < iV — 1 the functions Qmix, 9) as 
the final intersection point of the line Rh^,v^{x) + tSf^^ „ {&) with dB for t < 0. 
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Clearly, q„i is smooth on the set of (x, 6) for which the integrand doesn't vanish, 
since OB is smooth (see Figure [4]) . We then have a local regularized version of the 
weight function Wj given by 

Wreg j {x , 6) — cxp ( — / (j{x — t9, 9) dr ) 
(39) • ri exp (- / <j{qra{x, 9) + tS^I^^ (^), m dr) , 

m— 1 ^ 

for < j < A^— 1. By construction, Wregjix, 9) coincides with Wj{x^ 9) or Wj-i{x, 9) 
on the broken rays we've isolated, with the index depending on whether the broken 
ray intersects c. This allows us to rewrite the ballistic part of the normal operator 
as 

ballistic fix) 

= E / / H^{^'-3{x,0),9-j{x,9))wj{x,9)[wjf]{x + T^{x,9)9 + t9,9)dtd9 



§1 JR+ 



[ [ \a\^{qj{x,9),S^\^,{9))Wreg.j{x,9)[Wreg.jfKx + T^{x,9)9 + t9,9)dtd9. 
~^ JSi JR+ ' ' 



From here, we can make the change of variables s = < + r_ (x, 0). Then from Lemma 
2 of |FSU08| we have a classical pseudodifferential operator of order —1. 

Now we need to deal with N^^E,a.refiect- However, the approach is similar. Recall 
we have 



N N » » 

Na.E.a.reflectfix) ^ E / / ^ ) I^P] (^^-ji (a^, 6'), (x, 6'))Wj, (x, ( 



\w 



02 



/](^i2-ii {x,9)+ t9.j^^.j^ {x,9), 9j^^j^ {x,9)) dtdO 



N-l N-1 



I I [Mo^W^-n{x,e),9^,,{x,9))w,,{x,9) 



Jl=0 j2=0j2#jl ' 

(40) • K- J](z,,_,, (x, 9) + (x, 9), (x, 9)) dtd9. 

Again suppose that the broken ray ^xafio intersects a corner c at the mth reflection, 
where 1 < m < A^ — 1. We claim that in the first summation above, the terms with 
j2 = ™ vanish. To see this, note that if (x,0) is such that z-j-^{x,9)^9-j-^{x,9) is 
in the support of (1 — l[/)|ap, then reflecting j2 = m times from that point will 
put us in the small "bow-tie" region near c on which the corresponding line integral 
vanishes. So we rewrite the flrst term in (40) as 

N N „ „ 

E E / / [{l-'^u)\aWz^,Ax,e),9,,,{x,9))w,,{x,9) 

(41) • [wjj]{z,,^,,{x,9)+t9,.^_^^{x,9),9,^_^^{x,9))dtd9 



The important thing is that in (41 1 the reflection signature of Zj^^j^ relative to 
for each pair (ji,j2) is the same wherever a^(l — are supported. 
As in S4.2 for each pair of indices (ji,j2): we consider the reflection signature 
when going from {x,9) to {zj^_j^,9j^-j^). This is given by (hj^ ~hj^,Vj^ ~Vj^). So 
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as in ( 32 ) we define 



a _ yjij2 ~ ^ 



\y 31,32 -^1 

and we can rewrite 



.E .a,re fleet 
AT-l Af-1 

= E E 

Vji ,32 



[il-lu)\a\'] [qn 



yji ,32 



X, 



\y 31,32 



N-1 N-1 

+ E E 

jl=0 j2=0j2#jl 



[lu\a\ 



^reg,j2 



yji, 
lyji. 



y 31,32 ^1 

dy 



o-l 
' 1^31 ."ii 



V31 ,32 



|yjlj2 



(I31 



0-1 



In 
l%i 



Wh,. 



■^reg,jAy^v) f{y) 



dy 



N~l N-1 

E E 

31=0 32=0,32^31 



yji, 



(42) 



X, 



yji ,32 



\y 31,32 



Wreg,j2iyiV) f{y) 



dy 



\y3i,32 



Finally, we are ready to prove the result about the decomposition of Ng-^E.a- 

Proof of Proposition^ Recall that we have already shown in the discussion, see 
([22]), that 

,E ,oc, ballistic ,E,a^reflect- 

where N„^e, a, ballistic is a classical pseudo differential operator of order —1 that is 
elliptic on M (more generally on the subset of the cotangent space given by M'). 
So we have a parametrix Q of Na-,E.a,baiiistic of order 1 that is elliptic on A4. 

For the reflected part, we have a sum of terms of the form ([32]) in the case of 
a neighborhood of broken rays that don't have any corner points (except possibly 
as a starting point), or Na-^E,a,re fleet has the form of (42 1, both of which have C§° 
Schwartz kernels. Applying Q to both sides of ([22|) yields the result. □ 



Proposition 3. Under the conditions of Theorem^ without assuming that Ia-,E,, 
is injective, 

(a) one has the a priori estimate 

WfWmK) < C\\N^,E,af\\min) + Cs\\f\\H~Hn). Vs; 

(b) Ker/o-.E.a is finite dimensional and included in C°°{K). 
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Proof. The proof is essentially the same as that for ( |FSU08j . Proposition 3). To 
prove part (a), by Proposition [2] we have that 

,E ,a.re fleet f + Slf f + S2f. 

S2 has Schwarz kernel and by elliptic regularity we have 

ll/l|L2(/f) - Cs\\f\\H- = {n) < \\QNcr^E,af\\L^{K) < C*! I /I I ffi (Q) 

as desired. 

For part (b), if / G Kei Ia,E,a, then (Id + 52)/ = 0, and also 6*2 is a compact 
operator on L'^{K) with smooth kernel. By the Fredholm theorem, / lies in a finite 
dimensional subspace of C°°(f2). □ 

5. Reducing the Smoothness Requirement on a, a 

Given a choice of smooth cutoff a and a smooth a vanishing near dQ and such 
that Ia,E,a is injective, we would like to be able to perturb a and a in slightly 
and still have Na,E,a be injective. We do this according to the following modified 
version of ( |FSU08| . Proposition 4). 

Proposition 4. Assume that a, a are fixed and belong to . Let {a' ,a') he 0{5) 
close to [a, a) in . Then there exists a constant C > that depends on an a 
priori bound on the norm of (cr, a) such that 

(43) \\{Na',E,o.' - N„^E,c.)f\\HUn) < CSWfh^i^K)- 



Proof Recall that N„^E,a = N^,E,a,ballistic + N^,E,a,reflect, whcrC N^^E, a, ballistic IS 

a weakly singular integral operator with smooth kernel and N„ E,a,refiect is smooth- 
ing. By the discussion in the proof of Proposition 4 of |FSU08| . we have 

(44) ||iV. .E .a.ballistie Na',E,a'.ballistic\\H^{n) < Cd\\f\\L2(K)- 

Similarly, the Schwartz kernels of N„^E.a,refiect depend smoothly on a and a (see 
(32 1 and (42)), so the same estimate holds for this term. The completes the proof. 

□ 

Proof of Theorem^ First note that if Ia,E,a ■ L,^{K) — > L^(r_) is injective, then 
Na.E,a '■ Li^{K) H^{n) is also injective by an integration by parts. Next, since 
N a. E, a, ballistic IS elliptic ou , wc may apply a parametrix Q to get the equation 

QNa,E,af = f + Sif + QN^^E,a,reflectf- NotC that Si + QN„^E,a,reflect is COmpaCt. 

So by elliptic regularity we have the estimate 

\\f\\L^{K) < C'(ll^(T,£;,a/||//i(f2) + IKS'! + QN„ 

,E,a,reflect )/IIl2(0))- 

By Lemma 2 of |SU05j . the compactness of the second term on the right implies 
that 

< C"||A'"o-,£;,a/||ifi(0)- 

The other inequality is obvious. 

Finally, to prove part (b), we just use Proposition |4] to get 

II/IIl^(a') <C||iV,„,£^„o/|| 

< C\\N„^^ E,a'>f\\H^{n) ~ C\\{Ncro,E.aO - ^CT,£;,a)/||ffi(n) 

< C\\N,„,E,aof\\min) - CH\\f\\L2(K)- 

□ 
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6. Discussion 

In this work, we have utilized the machinery of |FSU08| for the X-ray transform 
with generic weights and apphed it to the given broken ray transform on the Eu- 
clidean square. The technique of path unfolding can be applied to any polyhedral 
domain, however the square is very nice in that reflections are well-defined even 
at corner points. The unfolding approach also placed significant restraints on the 
attenuation coefficient cr, which is why we could only reasonably tackle the case 
(7 = and then use the normal operator computation to extend the results to Cq 
perturbations. It is perhaps of interest to carry out similar computations for the 
broken ray transform on higher dimensional polyhedral domains, although it is not 
clear how one would be able to deal with corners other than to cut them out of the 
measurements altogether. One interesting problem to consider would be to try to 
get analogous microlocal results for the unit ball. However, one potential obstacle 
for strictly convex smooth domains is that a beam of collimated rays, upon reflec- 
tion becomes uncoUimated. This makes it difficult, or perhaps impossible, to use 
smooth cutoff functions of the measurements to isolate individual beams of broken 
rays from which one can obtain a nice integral form of the normal operator. Re- 
gardless, it is the author's hope that more progress can be made on this particular 
problem and its variants in the future. 
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